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ABSTRACT 

Most genomes, including yeast Saccharomyces 
cerevisiae, are pervasively transcribed producing 
numerous non-coding RNAs, many of which are 
unstable and eliminated by nuclear or cytoplasmic 
surveillance pathways. We previously showed that 
accumulation of PH084 antisense RNA (asRNA), in 
cells lacking the nuclear exosome component Rrp6, 
is paralleled by repression of sense transcription in 
a process dependent on the Hda1 histone 
deacetylase (HDAC) and the H3K4 histone methyl 
transferase Set1. Here we investigate this process 
genome-wide and measure the whole transcriptome 
of various histone modification mutants in a Arrp6 
strain using tiling arrays. We confirm widespread 
occurrence of potentially antisense-dependent 
gene regulation and identify three functionally 
distinct classes of genes that accumulate asRNAs 
in the absence of Rrp6. These classes differ in 
whether the genes are silenced by the asRNA and 
whether the silencing is HDACs and histone methyl 
transferase-dependent. Among the distinguishing 
features of asRNAs with regulatory potential, we 
identify weak early termination by Nrd1/Nab3/Sen1, 
extension of the asRNA into the open reading frame 
promoter and dependence of the silencing capacity 



on Set1 and the HDACs Hda1 and Rpd3 particularly 
at promoters undergoing extensive chromatin 
remodelling. Finally, depending on the efficiency of 
Nrd1/Nab3/Sen1 early termination, asRNA levels are 
modulated and their capability of silencing is 
changed. 

INTRODUCTION 

The development of high-density tiling arrays and 
large-scale RNA sequencing approaches revealed that all 
eukaryotic genomes are pervasively transcribed and syn- 
thesize a myriad of non-coding RNAs (ncRNAs). The 
yeast Saccharomyces cerevisiae produces a large number 
of long ncRNAs between 200 up to a few thousand bases 
in length, comprising both intergenic RNAs and tran- 
scripts antisense to coding open reading frames (ORFs). 
Antisense transcription can be regulated in a pre-initiation 
complex-dependent manner independently of divergent 
sense transcription of a downstream gene (1). Antisense 
transcripts (asRNAs) have been divided into stable 
unannotated transcripts (SUTs) detected in wild-type 
(WT) cells and cryptic unstable transcripts (CUTs), 
which are rapidly degraded and detectable only in 
mutants of the nuclear exosome (2,3). The degradation 
of CUTs by the 3—5' exonuclease Rrp6 is tightly linked 
to transcription termination by the Nrdl/Nab3/Senl 
(NNS) complex and polyadenylation by the TRAMP 
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complex (4-6). NNS is recruited at the 5' end of transcrip- 
tion units by RNA polymerase II (PolII) phosphorylated 
on serine 5 and induces termination primarily of 
<l-kb-long ncRNAs on recognition by Nrdl and Nab3 
of specific sequence motifs on the nascent transcripts 
(7-12). Recent global studies showed that Nrdl and 
Nab3 have a widespread localization and preferentially 
bind divergent and antisense ncRNAs playing a role in 
surveillance through early termination of transcripts that 
originate from bidirectional promoters or from the 3' end 
of ORFs (13-15). 

Growing evidence suggests that pervasive transcription, 
and in particular antisense transcription, contributes to 
the regulation of gene expression. Large-scale analyses 
of sense/antisense pairs, primarily considering SUTs, 
revealed anti-correlated expression patterns and anti- 
sense-dependent fine-tuning of sense expression (16). 
Gene-specific studies indicate that ncRNA transcription 
can modulate gene expression by mechanisms of transcrip- 
tion interference, nucleosome remodelling or via changes 
in histone modifications (17). Examples include the 
upstream ncRNA SRG1, which interferes with transcrip- 
tion of the downstream SER3 gene by enhancing nucleo- 
some assembly over the promoter (18-20), or the RME2 
asRNA, proposed to interfere with transcription elong- 
ation of the IME4 gene (21,22). H3K4 dimethylation de- 
posited by the H3K4 histone methyl transferase (HMT) 
Setl during antisense or upstream non-coding transcrip- 
tion has also been involved in gene repression by signalling 
the recruitment of the Rpd3 or Set3 histone deacetylase 
(HDAC) complexes (23-27), and genes repressed by Setl 
are enriched in antisense-producing ORFs (28). 

PH084 is another gene regulated by an antisense 
transcript. Our earlier studies indicated that the accumu- 
lation of PH084 asRNA in the absence of the exosome 
component Rrp6 is paralleled by repression of sense tran- 
scription through the recruitment of the Hdal/2/3 HDAC 
complex and targeted histone deacetylation at the PH084 
promoter and 5' end (29). We also showed that Setl 
promotes antisense production and contributes to 
PHOS4 gene repression (30). More recent single 
molecule fluorescence in situ hybridization (FISH) 
studies from our laboratory indicate that the expression 
of PH084 sense and asRNAs is anti-correlated in individ- 
ual cells. We also provided evidence that PH084 asRNA 
early termination by Nrdl is reduced in Arrp6 and that 
asRNAs escaping this early termination do not accumu- 
late in the nucleus but are rapidly exported. The data 
suggest that PH084 repression depends, at least in part, 
on constant low-frequency antisense transcription through 
the promoter, which is increased in Arrp6, reinforcing the 
tight on-off switch of this highly regulated gene (31). 

Although regulation by antisense transcription appears 
as a common mechanism (16), its molecular basis has not 
been examined globally. In this study, we aimed at 
characterizing the mechanistic diversity of antisense- 
mediated repression. To do so, we investigated the role 
of chromatin modifications in gene repression following 
asRNA accumulation in Arrp6. We used high-density 
tiling arrays to examine the contribution of the HMT 
Setl or the HDACs Hdal and Rpd3 in asRNA 



production and silencing on a genome-wide level. Our 
data identify three classes of antisense transcripts: two 
that are associated with gene repression and differ on 
whether the silencing mechanism involves histone-modify- 
ing activities and one class that does not affect sense tran- 
scription. We show that the repressive effect of antisense 
transcription on sense expression is linked to the efficiency 
of asRNA early termination by NNS. 

MATERIALS AND METHODS 

The yeast strains and sequences of all the primers used in 
this study are listed in Supplementary Table SI. 

Sample preparations and analyses 
RNA preparation and tiling arrays 

Total RNA was treated with RNase-free DNasel using 
Turbo DNA-free kit (Ambion). For first-strand comple- 
mentary DN A (cDN A) synthesis, 20 ug of total RNA was 
mixed with 1.72 ug of random hexamers, 0.034 ug of 
oligo(dT) primer and incubated at 70°C for lOmin 
followed by lOmin at 25°C and then transferred on ice. 
The synthesis included 2000U of Superscript II Reverse 
Transcriptase, 50 mM Tris-HCl, 75 mM KC1, 3mM 
MgC12, 0.01 M dithiothreitol, dNTP + dUTP mix 
(0.5 mM for dCTP, dATP and dGTP; 0.4 mM for 
deoxythymidine triphosphate and 0.1 mM for dUTP, 
Invitrogen), 20 ug/ml actinomycin D in a total volume of 
105 ul. The reaction was carried out in 0.2-ml tubes in a 
thermal cycler with the following thermal profile: 25°C for 
lOmin, 37°C for 30min, 42°C for 30min followed by 
lOmin at 70° for heat inactivation and 4°C on hold. 
Samples were then subjected to RNase treatment of 
20min at 37°C (30U of RNase H, Epicentre, 60U of 
RNase Cocktail, Ambion). First-strand cDNA was 
purified using the MinElute PCR purification kit 
(Qiagen), and 5 ug was fragmented and labelled using 
the GeneChip WT Terminal labelling kit (Affymetrix) ac- 
cording to the manufacturer's protocol. The labelled 
cDNA samples were denatured in a volume of 300 ul con- 
taining 50 pM control oligonucleotide B2 (Affymetrix) 
and Hybridization mix (GeneChip Hybridization, Wash 
and Stain kit, Affymetrix) of which 250 ul was hybridized 
per array (S. cerevisiae yeast tiling array, Affymetrix, PN 
520055). Hybridizations were carried out at 45°C for 16 h 
with 60rpm rotation. The staining was carried out using 
the GeneChip Hybridization, Wash and Stain kit with 
fiuidics protocol FS450_0001 in an Affymetrix Fluidics 
station. 

RNA analyses by real-time quantitative polymerase chain 
reaction 

Total RNA was prepared as described (29). For real-time 
quantitative polymerase chain reaction (RT-qPCR) quan- 
tifications, 1 ug of cDNAs of sense or antisense RNAs 
(asRNAs) was generated by Superscript II reverse tran- 
scriptase (Invitrogen) from total RNAs using strand- 
specific DNA primers. cDNAs were quantified by 
RT-qPCR (BioRad). The same amplicon was used to 
quantify sense and antisense cDNA. 
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Chromatin immunoprecipitation 

Chromatin immunoprecipitations (ChlPs) were performed 
essentially as described previously (29). Yeast strains 
were grown to OD600 = 0.8 either in yeast extract- 
peptone-dextrose (YEPD) growth media at 25° C and 
cross-linked for lOmin by the addition of formaldehyde 
to a final concentration of 1.2%. Cross-linked and 
sonicated chromatin extracts from 1.5 mg of Bradford- 
quantified proteins were immunoprecipitated overnight 
in the presence of protein G Sepharose (Amersham 
Pharmacia) with antibodies against H3K4me3 (Abeam 
8580) or H3 (Abeam 1791). All immunoprecipitations 
were repeated at least three times with different chromatin 
extracts. Immunoprecipitated DNA was purified and 
quantified by RT-qPCR with primers listed in 
Supplementary Table SI and expressed as the percentage 
of input DNA normalized to H3. Error bars correspond 
to standard deviations. 

Half-life measurements 

Cells were grown to an OD 600 = 0.8 in YEPD medium. At 
T = 0, 100(ig/ml 1,10-phenanthroline (Sigma) was added 
to the culture as described (32), and samples were collected 
at different time points and analysed by RT-qPCR. 
Half-lives were calculated by the equation t lj/2 — 0.693/k, 
where k is the rate constant for messenger RNA (mRNA) 
decay. Values of each time point are normalized for 
internal variations with SCR1 RNA, a control that is 
still stable at the 30-min time point. 

Anchor away strains 

The anchor away strains (Rrp6-AA and Nrdl-AA) were 
constructed by tagging the endogenous gene of interest at 
the carboxyl terminus by PCR amplification of the 
pFA6a-FRB-KanMX6 cassette (see Supplementary 
Table SI for primer sequences) and transformation of 
the PCR product into the parental rapamycin-insensitive 
strain HHY168 (33). Cells were grown to an OD 600 = 0.8 
in YEPD medium, treated with 1 ug/ml rapamycin 
(LC laboratories) and collected at different time points. 
RNA was extracted and analysed by RT-qPCR (or 
northern blotting). Values of each time point are 
normalized to ACT1 mRNA, a control not affected by 
rapamycin treatment. 

Bioinformatics analyses 
Gene classification 

Our goal was to classify genes according to whether they 
are regulated by antisense transcription and/or chromatin 
modification enzymes. To capture the effect of antisense 
transcription on the sense expression simultaneously for 
all histone modification mutants, we decided to use an 
unsupervised clustering approach in which we classify 
genes based on the combined effect of all mutants on 
the genes sense and antisense expression levels. For this, 
we calculated differential expression values in Arrp6 versus 
WT, AsetlArrp6 versus Arrp6, Ahda2Arrp6 versus Arrp6 
and Arpd3Arrp6 versus Arrp6 for the sense transcript from 
— 100 to +750 bp of the transcription start site (TSS) and 
for the antisense transcript —100 to +750 bp of the tran- 
scription termination site (TTS) on the antisense strand in 



bins of 20 bp. We then concatenated these values for each 
gene and performed pam (partitioning around medoids) 
clustering (34) using 15 clusters (Supplementary Figure 
SI). To ensure the robustness of the clustering, we per- 
formed 120 additional clusterings with 10-15 clusters 
(20 for each) randomly selecting 80% of all genes and 
measured how often each gene is paired with any other 
gene of the same cluster. Looking at the distributions, we 
selected a cut-off of 60% to identify the core set of genes 
for each cluster (Supplementary Figure SI). We then in- 
spected the median profiles of the 1 5 clusters and grouped 
them into five classes based on their expression profiles 
across the different conditions as follows: Class (I): 
cluster 1; Class (II): cluster 2; Class (III): clusters 6, 8 
and 12; Class (IV): cluster 5; and unclassified control 
set: clusters 3, 4, 7, 9, 10, 11, 13, 14 and 15 
(Supplementary Figure SI). The classes are described in 
Figure 2 and the corresponding genes listed in 
Supplementary Table S2. We focused here on characteris- 
tics that were interesting in the context of this study. 
Therefore, we did not consider clusters that showed 
elevated expression levels of antisense transcripts only in 
some of the double mutants but not in the Arrp6 versus 
WT separately (e.g. elevated antisense in Arpd3Arrp6 
versus Arrp6 in clusters 9 and 14) (see Supplementary 
Figure SI). Tiling array heatmaps can be viewed at 
http : / /steinmetzlab .embl . de//cgi-bin /viewStutzLab Array . 
pl?showSamples = stutzLabArray&type = heatmap&gene 
= rrp6. 

The complete tiling array data have been submitted to 
Array Express (at http://www.ebi.ac.uk/arrayexpress/, 
accession number E-MTAB-1316). 

Expression levels of antisense RNA 

To calculate the mean expression profile of the antisense 
transcripts for the different classes, we obtained the 
expression level in Arrp6 on the antisense strand and 
took the median across the genes in a given class for 
each bin of 10 bp. We then aligned this matrix once to 
the ORFs TSS and once to its TTS and calculated a 
moving average with a window size of 150 bp and a step 
size of 10 bp for each class across all genomic positions. 

Promoter structure, Nrdl/Nah3 motifs and Nrdl binding 

To analyse the promoter structure, we used the definition 
of TATA boxes from Basehoar (35) and the open-closed 
promoter configuration from Zaugg (36). TF binding sites 
were taken from Maclsaac (37). 

To analyse the occurrence of Nab3 and Nrdl motifs in 
antisense transcripts, we used the consensus sequences 
TCTT and GTA[AG] for Nab3 and Nrdl, respectively 

(12) . For each gene, we then counted the number of 
motifs within 400-0 bp upstream of the genes TTS on 
the antisense strand. For analysing Nrdl binding events, 
we obtained PAR-CLiP data from the Corden laboratory 

(13) . We summed the binding data within 400-0 bp 
upstream of the genes TTS on the antisense strand and 
normalized it to the WT expression level of this region. 
Using expression level in Arrp6 for normalization did not 
change the results. 
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Identification of the whole population of ncRNAs 

To identify specific ncRNA transcripts affected by Setl 
(Supplementary Figure S5) we performed an automatic 
segmentation of the tiling array data in the Arrp6 
mutant (38). We then overlapped the expressed segments 
with annotated transcripts identified by Xu et al. (3). 
Transcripts that did not overlap with the Xu et al. data 
and were at least 200 bases long were annotated as 
ncRNAs (Supplementary Table S3). We redefined CUTs 
and SUTs such that CUTs are at least 2-fold upregulated 
in Arrp6 versus WT (Supplementary Table S3 and 
Supplementary Figures S7A and S7B). We found more 
ncRNAs than Xu et al. (2009), likely because we used 
total RNA instead of poly(A)-enriched RNA and per- 
formed the experiment in a different strain background 
(W303 instead of S96/S288c). This annotation has been 
used for the analyses in Supplementary Figures S5 and S7. 

Comparison of H3K4me3 levels in A setl -silenced versus 
not affected genes 

We grouped the ORF transcripts as well as the ncRNA 
transcripts into four quartiles, each based on their expres- 
sion level in Arrp6. Transcripts were then separated into 
'Asetl-silenced' and 'not affected' based on their differen- 
tial expression value in AsetlArrp6 versus Arrp6, with 
thresholds of 2-fold over- or under-expressed 
(Supplementary Figure S5D). H3K4me3 in their 
promoter region (—100 to +300 bp of the TSS) was 
measured using data from (39). 

RESULTS 

A clustering approach identifies different classes of genes 
potentially regulated by asRNAs 

To identify genes regulated by antisense ncRNAs 
accumulating in the absence of Rrp6 and to characterize 
the mechanisms by which antisense transcription affects 
gene expression, we performed tiling array expression 
profiling of strains lacking Rrp6 alone or in combination 
with disruptions of histone modification enzymes previ- 
ously implicated in antisense-mediated regulation of 
PH084: the H3K4 methyl transferase Setl and the 
Hdal/2/3 HDAC complex subunit Hda2. In addition, 
we selected the HDAC Rpd3, which has been involved 
in antisense-mediated regulation of other genes (23,24). 
Total RNA prepared from single and double mutant 
strains in triplicate was hybridized to strand-specific 
high-resolution tiling arrays covering the 5. cerevisiae 
genome as described (3,38). 

We first validated the mutant phenotypes by examining 
the behaviour of our model gene PH084 by quantitative 
RT-PCR (Figure 1). We confirmed PH084 repression in 
response to antisense accumulation in Arrp6 and its de- 
repression on loss of Hda2, a component of the Hdal 
complex. In addition, we found a similar de-repression 
of the gene on loss of the HDAC Rpd3 in Arrp6, suggest- 
ing that the Hdal and Rpd3 HDACs play partially 
overlapping roles in antisense-mediated silencing. 
Notably, the increase of PH084 sense RNA levels in 
Ahda2 or Arpd3 single mutants indicates a repressive 
effect of these HDACs in WT cells. This repression 
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Figure 1. Analysis of PH084 sense and antisense transcripts in differ- 
ent chromatin-modifying mutants. PH084 sense and asRNA levels 
were quantified by RT-qPCR in the indicated strains as described 
(29,30). Sense and antisense values were normalized to ACT1 and ex- 
pressed as fold change over WT. Error bars reflect standard deviations 
of an average obtained from three independent experiments. 

could be linked to the low levels of PH084 antisense ex- 
pression in WT cells (29,31). Finally, we confirmed that 
loss of the HMT Setl reduces antisense levels and PH084 
repression in Arrp6, supporting a role of Setl in silencing 
by acting upstream of asRNA, either by enhancing its 
production or stability. 

To study the potential silencing mechanisms mediated 
by Rrp6-sensitive asRNAs genome-wide, we used a clas- 
sification approach that integrates a gene's response to 
various chromatin mutants simultaneously for its sense 
and antisense transcripts. Briefly, we used strand-specific 
tiling array data to compute differential expression in 
sense and antisense directions for Arrp6 versus WT, 
Arrp6Asetl versus Arrp6, Arrp6Ahda2 versus Arrp6 and 
Arrp6Arpd3 versus Arrp6, and classified genes based on 
the behaviour of their sense and antisense transcripts in 
those conditions. This integrative approach allowed us to 
detect asRNAs whose differential expression falls below 
the detection threshold in individual conditions, as long as 
they show differential expression above 0 in any of the 
other mutant strains. The differential expression of sense 
mRNAs was calculated between —100 and +750 bp 
around the ORF TSS, whereas the differential expression 
of asRNAs was calculated on the opposite strand between 
-100 and +750 bp around the ORF TTS. We then applied 
unsupervised clustering using the concatenated differential 
expression values as features to group genes with similar 
sense/antisense behaviours across the different strains 
(see 'Materials and Methods' section, Supplementary 
Figure SI). A total of 1715 genes in convergent configur- 
ation that overlap each other were excluded from the 
analysis. To further corroborate the clustering, we per- 
formed a robustness analysis and selected only the core 
set of genes that co-segregated in multiple randomized 
rounds of clustering, strengthening the assessment of 
their distinctive features (Supplementary Figure S2A). 

These clusters were then grouped into four classes and 
a group of unclassified genes based on a combination 
of three criteria: first whether they exhibit asRNA 
accumulation in Arrp6, second whether the gene (ORF) 
itself is repressed in Arrp6 and finally whether the 
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Figure 2. Gene classification based on the silencing mechanism associated with Rrp6-sensitive asRNAs. Genes were clustered into 15 groups based 
on differential expression of the sense and antisense strands in four conditions: Arrp6 versus WT, Asetl Arrp6 versus Arrp6, Ahda2Arrp6 versus Arrp6 
and Arpd3Arrp6 versus Arrp6 (see Materials and Methods' section, Supplementary Figures SI and S2). The clusters were grouped into four Classes, 
(I), (II), (III) and (IV), based on the similarity of their median profiles for the different mutants and an unclassified control set (Unclassified). The 
number of genes in each class is given at the top. The distinguishing features for each class are highlighted in yellow and represented by little 
schematics where red indicates up- and blue indicates downregulation in a given mutant. Genes where silencing involves Hdal, Rpd3 and Setl are 
described as 'HDAC/HMT-dependent silencing'. 
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zl?77;<5-dependent repression of the gene is relieved on the 
loss of histone modification enzymes (Figure 2 and 
Supplementary Table S2). A class was defined as 
accumulating or silencing transcripts in a given strain if 
the median differential expression of the class was above 
or below 0, respectively. Class (I) contains genes that are 
repressed and accumulate asRNA in Arrp6, and de-re- 
pressed on the loss of Hda2, Rpd3 or the HMT Setl. 
The asRNAs in this class are essentially not influenced 
by these HDACs and modestly affected by the HMT 
Setl. Therefore, this class represents genes potentially 
regulated by an antisense-mediated mechanism similar to 
PHOS4. Genes in Class (II) also accumulate asRNA and 
are repressed in Arrp6. However, unlike Class (I) genes, 
the loss of Setl, Hda2 and Rpd3 only marginally affects 
ORF expression, indicating that potential antisense- 
mediated silencing is largely independent of these 
HDACs and the Setl HMT. Class (III) genes also show 
increased antisense accumulation in Arrp6 but no effect on 
ORF transcription. Notably, genes in both Classes (II) 
and (III) show substantial decrease of antisense transcripts 
in the absence of Setl, rendering these classes of genes the 
most dependent on Setl for antisense production in Arrp6. 
Class (IV) contains genes that are repressed by Setl and 
the two HDACs. These genes are related to Class (I) (see 
below); however, they segregate as a separate group, as 
they show either no clear asRNA signal within 750 bp 
upstream of the TTS (in some cases, asRNAs start close 
to the 5' end of the gene, e.g. PDR15), or asRNA levels 
that only modestly change in Arrp6 (similar to SUTs) and 
are not paralleled by changes in sense RNA levels (e.g. 
AT02), or sense transcription is mostly repressed under 
these experimental conditions (e.g. CRF1, YTP1). Finally, 
the remaining genes were labelled as 'unclassified 1 and 
used as a control set throughout the study. Notably, 
single mutant versus WT comparisons of the different 
gene classes, in particular Class (I), do not show substan- 
tial increases in ORF expression, supporting the view that 
the repressive effects of Setl, Hda2 and Rpd3 on these 
genes in Arrp6 are linked to the presence of asRNA 
(Supplementary Figure S2B). Representative expression 
patterns of genes from Classes (I) to (III) are shown in 
Supplementary Figure S3. 

In summary, we found ~100 genes [Class (I) and (II)] 
that are silenced concomitant to asRNA accumulation in 
Arrp6. While the proposed antisense-mediated repression 
of Class (I) (28 genes) depends on the Hdal/Rpd3 HDACs 
and the Setl HMT, silencing of Class (II) (69 genes) 
appears to involve a different mechanism possibly based 
on transcription interference (see below). However, we 
cannot exclude that the slight decrease in a fraction of 
Class (I) and (II) asRNA levels observed in Arrp6 in the 
absence of Setl, Rpd3 or Hda2 is due to increased sense 
transcription in these mutants or that the chromatin modi- 
fiers independently affect the two transcription units. 
Class (III) contains 469 genes that accumulate asRNAs 
in Arrp6 with no effect on gene expression. In addition, 
Class (IV) contains 76 genes that are repressed by the two 
HDACs and the HMT already in WT cells. Many of them 
produce low, mostly Rrp6-insensitive levels of asRNAs, 
making them distinct from, yet related to, Class (I). 



Finally, 2024 genes that have essentially no asRNA and 
are affected in none of the mutants remained unclassified. 

Functional and non-functional asRNAs show distinct 
features 

To gain further insights into antisense-mediated regulatory 
mechanisms, we investigated general functional features of 
the sense/antisense pairs in each class. We found that genes 
regulated by Rrp6-sensitive asRNA [Classes (I) and (II)] 
have significantly lower basal expression levels in WT 
compared with those that are not [Class (III)] 
(Figure 3A). This is in agreement with previous studies re- 
porting that genes affected by antisense SUTs tend to be 
expressed at medium to low levels (16). Interestingly, Class 
(IV) genes, which are de-repressed in the HDAC/HMT 
knockouts, are lowly expressed, consistent with a generally 
repressed state in current growth conditions. Accordingly, 
promoters of genes in this class are enriched for binding 
sites of repressive transcription factors (TFs), such as 
glucose repression (MIG1) and meiosis (UME6, IME1), 
but also of several activators (ADR1, MSN2, MSN4, 
SPT15/TBP and STP1), indicating extensive regulation 
(Supplementary Table S4 and see below). 

To further define the features of antisense transcripts that 
affect the corresponding sense gene (termed 'functional 
asRNAs'), we compared the extent to which they overlap 
with the sense gene as well as their expression levels with 
asRNAs that have no effect on sense transcription (termed 
'non-functional asRNAs'). For this, we calculated the 
average expression in Arrp6 along the antisense strand sep- 
arately for each class. We aligned the expression values once 
to the TTS of the corresponding sense ORF and once to its 
TSS (Figure 3B). Interestingly, despite their relatively low 
expression level, asRNAs in Class (I) extend across the 
whole length of the gene and even beyond the TSS. In 
contrast, the expression of Class (III) antisense transcripts, 
which have no effect on sense expression, shows a strong 
peak close to the TTS but then drops rapidly, indicating 
that these transcripts do not extend much into the genes. 
Class (II) antisense transcripts exhibit an intermediate 
pattern, as they show an initial peak like in Class (III), 
but decrease more slowly and appear to extend beyond 
the TSS, consistent with an effect of these asRNAs on 
ORF expression (Figure 3B). The asRNA profile for 
Class (IV) genes shows a similar pattern to Class (I) genes 
but at a much lower level, while no asRNA was detected 
over the unclassified control set. Comparing the length of 
the genes/ORFs across the different classes revealed no 
specific bias that could explain the difference in the extent 
of antisense overlap (data not shown). 

To follow up on the finding that Class (I) and 
(II) asRNAs tend to extend into the promoters of their 
counterpart genes, we examined whether there are any dif- 
ferences in promoter structures, such as TATA boxes or 
promoter-nucleosome architecture across the five classes. 
TATA box genes mainly comprise inducible or stress- 
response genes (35). Closed promoters, characterized by a 
well-positioned +1 and a fuzzy —1 nucleosome that often 
occlude important TF binding sites, have been associated 
with a complex regulatory structure and high variation of 
expression levels across a population of cells (36). We found 
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that Class (I) is the only class enriched for closed promoters 
and both Class (I) and Class (IV) are enriched for TATA- 
box promoters (Figure 3C). Consistently, they also exhibit 
significantly higher expression noise (40) compared with the 
other gene classes (Figure 3D). 

In summary, our data suggest that the ability of asRNAs 
to participate in gene repression resides in the extent to 
which the ncRNA transcription progresses along the gene 
as well as in the corresponding promoter structure. The 
enrichment of Class (I) genes for closed promoter config- 
urations suggests that silencing could involve antisense- 
induced nucleosome repositioning in the sense promoter, 
whereas another mechanism, less dependent on chromatin 
remodelling, may be operating at Class (II) genes. 

Setl is involved in asRNA production 

Our clustering analysis revealed that a large fraction of the 
asRNAs that accumulate in Arrp6, especially those of 
Class (II) and Class (III), tend to be downregulated on 



deletion of the HMT Setl (Figure 2, AsetlArrp6 versus 
Arrp6). Setl is the catalytic subunit of the COMPASS 
complex devoted to the differential methylation of lysine 
4 on H3 in S. cerevisiae. H3K4 trimethylation by Setl is 
enriched at the TSS of transcribed ORFs, while di- and 
monomethylated H3K4 peak in the middle and 3' end of 
the gene, respectively (39,41). 

We next asked whether Setl activity is also associated 
with transcription of asRNA. For this, we compared the 
genome-wide distribution of H3K4me3 marks at the TTS 
of genes in different classes using published ChlP-chip 
data (39). We found that the H3K4me3 levels between 0 
and — 300 nt from the TTS of ORFs that have a clear 
asRNA [Classes (I), (II) and (III)] are significantly 
higher compared with those at ORFs with low or no anti- 
sense [Class (IV) and unclassified control set] (Figure 4A), 
which is in agreement with recently published observa- 
tions (28). ChIP analysis of H3K4me3 profiles on four 
specific genes (39), i.e. PH05 and GYP5 (Class I), RRI1 
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(Class II) and RAD4 (Class III), confirmed increased levels 
of this modification at their 3' ends (Figure 4B). To further 
support the hypothesis that the H3K4me3 peak at the 3' 
end of antisense-producing genes reflects antisense tran- 
scription, we examined H3K4 trimethylation over the 



Class (I) genes GYP 5 and PH05 in WT and Apho4 
strains (Figure 4C). Pho4 is a TF dispensable for GYP 5 
transcription but essential for the expression of several 
PHO genes, including PH05 and PH084. In Apho4, 
H3K4 trimethylation peaks were detected both at the 5' 
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and 3' end of the Pho4-independent G YP5 gene, while for 
PH05 only a single peak of H3K4me3 was detected at the 
3' end, consistent with the absence of sense but mainten- 
ance of antisense transcription (see Supplementary Figure 
S4A for PH05 and GYP5 sense and asRNA levels in these 
strains). As a control, we performed ChIP on the anti- 
sense-deficient ACT1 gene, where, as expected, we 
detected a single peak at the 5' end. Notably, the level of 
H3K4me3 at the 3' end of both GYP5 and PH05 was 
significantly enhanced on loss of Rrp6 in both WT and 
Apho4 backgrounds, indicating an increase of H3K4 
trimethylation by Setl in the absence of the exosome 
(Figure 4C). 

To link the effect of Setl on asRNA production to 
H3K4 methylation, we compared the levels of different 
ncRNAs in Asetl and the H3K4A mutant strain. In a 
Arrp6 background, the expression levels of asRNAs be- 
longing to Class (I) (PH05 and GYP5), Class (II) (RRI1) 
or Class (III) (RAD4) were similarly reduced in combin- 
ation with Asetl or H3K4A compared with Arrp6 alone 
(Figure 4D and Supplementary Figure S4B), indicating 
that the positive effect of Setl on ncRNA production 
depends on the methylation of H3 on lysine 4. asRNA 
decay measurements in Asetl and K3K4A mutant 
strains showed similar half-lives compared with WT, sup- 
porting the view that H3K4 methylation by Setl contrib- 
utes to asRNA production rather than stability 
(Supplementary Figure S4C). Replacing Asetl by a Setl 
catalytic mutant (setlca, G951S) (42) in our tiling array 
analyses had slightly weaker, yet similar, effects on 
asRNA expression (Spearman correlation if = 0.71 and 
P <2.2e-16) as well as on ORFs (Spearman correlation 
R = 0.63 and P<2.2e-16) (Supplementary Figures S4D 
and see below). These observations together indicate 
that in addition to H3K4 methylation, the integrity of 
the Setl complex may also contribute to the production 
of asRNAs. 

We next sought to investigate the genome-wide contri- 
bution of Setl to ncRNA production. Previous studies on 
global transcriptional effects of H3K4 methylation 
focused on the analysis of ORF mRNA expression, and 
most identified little contribution of Setl to gene expres- 
sion (28,43^16). We first compared ORF and antisense 
differential expression in AsetlArrp6 versus Arrp6 versus 
Arrp6 versus WT. These analyses revealed a negative cor- 
relation for the asRNAs (—0.43), while no evident correl- 
ation was present for the ORFs (—0.096) (Figure S5A). 
This indicates that loss of Setl mainly affects the expres- 
sion of antisense transcripts that accumulate on loss of 
Rrp6. Additional analyses of the Asetl and setlca 
strains using the whole population of ncRNAs (see 
Materials and Methods' section for definition and 
Supplementary Table S3) confirmed that Setl preferen- 
tially affects CUTs versus SUTs, and is generally more 
important for the expression of ncRNAs than ORFs 
(Supplementary Figure S5B and C). 

To further understand the role of Setl, we studied 
the global relationship between H3K4me3 levels and 
Setl -dependent differential expression. For this, we 
compared H3K4me3 levels in WT (39) between ORF 
transcripts and ncRNAs that are downregulated in 



Asetl Arrp6 and those that are not affected by loss of 
Setl (Supplementary Figure S5D). Because H3K4me3 
correlates strongly with expression levels (39), we 
grouped genes and ncRNAs into four groups each 
based on their expression levels and analysed them sep- 
arately. We found that a large fraction of ncRNAs 
(>30%) and only a modest subset of genes (<10%) are 
affected by loss of Setl. The affected ncRNAs and ORFs 
show increased H3K4 trimethylation compared with non- 
affected transcripts, and this increase is most significant 
for lowly expressed ORFs or ncRNAs. These data 
confirm a specific role of Setl and H3K4me3 in 
promoting asRNA production. 

In summary, we showed that antisense production is 
associated with H3K4 methylation by Setl. Strikingly, 
the H3K4 trimethylation marks at the 3' end of asRNA- 
producing genes increase on loss of exosome activity 
(Figure 4C), suggesting that asRNA accumulation in 
Arrp6 might not be due to only stabilization but also 
increased transcription elongation, a process stimulated 
by Setl -dependent H3K4 methylation. These data 
extend our previous observations on PH084 (31) to a 
more general mechanism of many antisense-producing 
genes, and further support the view that loss of Rrp6 
facilitates escape from early termination. 

Early termination by the NNS complex differentially 
modulates AS accumulation and silencing capability 
in the different classes 

Rrp6-sensitive ncRNAs (or CUTs) are normally 
eliminated in a process that involves transcription 
termination by the NNS complex. To better understand 
the role of early termination in the different classes, we 
next investigated Nrdl and Nab3 binding across the func- 
tional [Classes (I) and (II)] and non-functional [Class (III)] 
antisense transcripts. Because the number of Nrdl and 
Nab3 binding motifs within the first few hundred base 
pairs correlates with early termination and degradation 
efficiency (5), we searched for Nrdl and Nab3 binding 
motifs on the antisense strand within 400 nt of the ORFs 
TTS. We then compared the number of motifs across 
functional and non-functional asRNAs. Interestingly, 
non-functional antisense transcripts from Class (III) 
genes contain significantly more Nrdl and Nab3 motifs 
than those that are functional [Classes (I) and (II) 
(Figure 5Ai)]. We also compared the Nrdl binding from 
PAR-CLiP data (13), and again found that non-functional 
asRNAs tend to be bound by Nrdl more often than 
functional asRNAs (Figure 5Aii). 

To address whether Class (III) genes may undergo more 
potent asRNA early termination by NNS than Class (I) 
genes, we examined the levels of a Class (I) (G YP5) and a 
Class (III) (RAD4) asRNA by northern blotting and 
RT-qPCR when depleting the essential Nrdl protein 
using the glucose-repressible GAL1 promoter in the 
presence or absence of Rrp6 (Figure 5B and Supplemen- 
tary Figure S6A). Both GYP 5- and RAD4-extended 
asRNA species strongly accumulate on depletion of 
Nrdl. Notably, the increase of GYP5 asRNAs is compar- 
able with that observed in Arrp6, while the levels of RAD4 
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asRNA are clearly higher than in Arrp6 (compare lanes 2 
and 3, GAL-NRDl in Glu and GAL-NRDl Arrp6 in Gal). 
These observations indicate that asRNAs with strong 
Nrdl/Nab3 termination signals are less sensitive to the 
absence of Rrp6. Northern blot analyses revealed no 
major differences in the size of asRNAs accumulating on 



Rrp6 or Nrdl depletion, indicating that similar termin- 
ation signals are used in the two conditions (Figure 5B). 
Our northern blot analyses were not able to detect the 
short early terminated and unstable asRNAs (CUTs) 
identified earlier in Arrp6 at the GYP 5 and RAD4 genes 
by 3'-long SAGE analysis (2), possibly because they are 
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too short and/or too few. Of note, the accumulation of 
GYP5 and RAD4 asRNAs detected on depletion of Nrdl 
was not significantly enhanced in the absence of Rrp6 
(lane 4), suggesting that these extended asRNAs are not 
subject to degradation by Rrp6, likely because they are 
exported into the cytoplasm (Castelnuovo,M., unpub- 
lished data). Importantly, while loss of Rrp6 only affects 
GYP5 sense expression, Nrdl depletion was accompanied 
by a reduction of both Class (I) {GYP 5) and Class (III) 
(RAD4) gene mRNA levels, indicating that abrogation of 
early termination can convert a non-functional into a 
functional asRNA (Figure 5B and Supplementary Figure 
S6A). This effect is unlikely due to direct modulation of 
GYP 5 and RAD4 mRNAs by Nrdl based both on PAR- 
CLiP analyses and the recent observation that NNS- 
mediated mRNA attenuation is a rare event (13,14). 

To more accurately investigate the kinetics of antisense 
accumulation and sense repression, we induced rapid de- 
pletion of Rrp6 or Nrdl from the nucleus using the 
rapamycin-inducible anchor away (AA) technique (33), 
and measured GYP5 and RAD4 sense and asRNA levels 
at various times after rapamycin treatment (Figure 6A and 
Supplementary Figure S6B). We confirm that loss of Rrp6 
and Nrdl results in a comparable increase in GYP5 
asRNA levels (5-fold versus 4-fold) accompanied by 
sense repression, while the accumulation of RAD4 
asRNA is weaker in Rrp6-AA compared with Nrdl-AA 
(8-fold versus 30-fold). Similar to the GAL-NRD1 deple- 
tion experiment (Figure 5B), RAD4 asRNA accumulation 
in the Nrdl-AA but not in the Rrp6-AA strain was 
accompanied by repression of sense transcription. 
Notably, during Nrdl nuclear depletion, both GYP5 and 
RAD4 asRNAs accumulate 20-30 min before sense repres- 
sion. In addition, there is an almost linear correlation 
between the amount of asRNA and the extent of silencing. 
We propose that sense repression depends on the 
frequency of elongated asRNAs escaping from an early 
termination control. 

asRNA-dependent repression of Class (I) genes such as 
PH084 or GYP 5 involves the HDACs Hda 1/2/3 and 
Rpd3 (Figures 1 and 2) and (29). To address whether 
the sense repression of Class (I) and Class (III) genes 
observed when rapidly depleting Nrdl requires histone 
deacetylation, we deleted Hda2 in the Nrdl-AA strain. 
Interestingly, repression of GYP5 and RAD4 sense tran- 
scription observed in the Nrdl-AA strain grown in 
rapamycin was completely rescued in the absence of 
Hda2 (Figure 6B). As previously proposed (29,31), the 
Hda 1/2/3 HDAC may trigger promoter deacetylation to 
consolidate repression initiated by antisense transcription 
through this region. These data suggest that asRNAs or 
antisense transcription may recruit histone-modifying 
enzymes to specific target sites as has already been 
proposed in higher eukaryotes (47). 



DISCUSSION 

Features of functional Rrp6-sensitive asRNA 

Our study addressed the genome-wide effects of antisense 
transcripts that accumulate in Arrp6, extending the single 



gene studies on PH084, a model gene for asRNA-mediated 
transcriptional silencing in Arrp6. We quantified the genes 
that become repressed through antisense transcription in a 
mechanism involving Setl and/or the HDACs Hdal and 
Rpd3, and identified major features of functional antisense 
transcripts (Figure 1). Among the 566 asRNAs that accu- 
mulate in Arrp6 [Classes (I), (II) and (III)], 97 [Classes (I) 
and (II)] are associated with the repression of sense gene 
expression (Figure 2). Earlier studies addressed the global 
effect of asRNAs that were classified as stable transcripts 
(SUTs) (16,48), which constitute a fraction of our 
Rrp6-sensitive asRNAs (Supplementary Figure S7A and 
B). Interestingly, our results suggest that functional 
asRNAs, whether previously studied SUTs or our Rrp6- 
sensitive transcripts, have similar features: both extend into 
the ORF promoter, and their repressive effects are targeting 
stress and other highly regulated genes (Figure 3) and 
(16,48). 

In addition, we found that the repressive effect of 
Rrp6-sensitive asRNAs inversely correlates with the 
efficiency of early transcription termination by the NNS 
complex. Our recent studies on PH084 provided evidence 
that NNS function is affected in Arrp6 (31). Based on the 
described physical interactions between NNS and the 
nuclear exosome (6), loss of Rrp6 may compromise early 
termination by weakening optimal Nrdl recruitment and 
binding to the transcript. Alternatively, loss of Rrp6 may 
indirectly affect early termination through accumulation 
of ncRNAs bound to Nrdl/Nab3, preventing the recycling 
of the complex. This Nrdl/Nab3 sequestration in Arrp6 
may primarily affect transcription attenuation of asRNAs 
with weak Nrdl/Nab3 binding sites. Accordingly, the fold 
increase in asRNA levels on depletion of Nrdl versus 
Rrp6 is weaker for Class (I) compared with Class (III) 
asRNAs (Figures 5 and 6A). Class (I) asRNAs, including 
PH084 antisense, but also Class (IV) asRNAs are sub- 
jected to weak early termination and can often be detected 
in WT cells (31) (data not shown). SUTs, which are less 
prone to early termination and degradation by Rrp6 
contain fewer Nrdl/Nab3 binding motifs than CUTs, 
and this feature is even more pronounced in ORFs 
(Supplementary Figure S7C). 

Our observations on asRNA accumulation following 
Nrdl depletion have recently been extended to the whole 
genome by the Cramer laboratory (14). This study 
identified 1526 Nrdl unterminated transcripts or NUTs 
on rapid nuclear depletion of Nrdl, which overlap with 
a large fraction of CUTs (68%) as well as SUTs (58%) but 
poorly with mRNAs. Thus, early termination by Nrdl 
ensures proper gene expression mainly by restricting 
divergent transcription from bidirectional promoters as 
well as asRNAs at most yeast ORFs. 

Role of Setl and Rrp6 in the production of asRNAs 

Despite the general occurrence of H3K4 methylation on 
active genes, we found, consistent with previous studies, 
that loss of Setl has only minor global effects on steady- 
state gene expression (28,43,44,46). However, we identified 
a role for Setl in promoting ncRNA production, a feature 
mostly revealed in Arrp6 (Supplementary Figure S5), when 
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Figure 6. (A) Anchor away of Rrp6 and Nrdl from the nucleus promotes asRNA production and gene repression. GYP5 and RAD4 sense (S) and 
antisense (AS) RNA levels were analysed by RT-qPCR in the anchor away strains Rrp6-AA (grey bars) and Nrdl-AA (black bars) grown in the 
presence of rapamycin (+) for the indicated times (minutes). The white bars correspond to S and AS RNA levels in the absence of rapamycin (— ). All 
values are expressed as fold change compared with the Rrp6-AA strain at time Omin and normalized to ACT1. (B) Hda2 contributes to Class (I) and 
(III) gene repression on anchor away of Nrdl. Northern blot analysis of GYP5 and RAD4 sense and antisense (AS) RNAs in the Nrdl-AA and 
Nrdl-AAzl/iJai strains grown either in the absence (— ) or presence (+) of rapamycin (RAP) for the indicated time in minutes. ACT1 was used as a 
loading control. The GYPS and RAD4 mRNA signals normalized to ACT1 were expressed as a percentage of the signal detected at time 0 in the 
same strain. 



early termination by Nrdl is compromised. Support for 
this role comes from a recent study, which proposed that 
while H3K4me2 is involved in sense repression, 
H3K4me3, deposited by Setl, contributes to asRNA 
production (28). Our results from half-life analyses in 
Asetl and H3K4 mutants further indicate that the 
positive effect of H3K4me3 on asRNA levels in Arrp6 is 
due to increased asRNA production rather than increased 



stability (Supplementary Figure S4C). Consistently, we 
showed that loss of Rrp6 alone results in increased 
asRNA-dependent H3K4 methylation over the PH084 
gene (31) and two other Class (I) genes PH05 and 
GYP 5 (Figure 4C). 

The preferential role of Setl in promoting Rrp6-sensi- 
tive asRNAs indicates a functional link between Setl and 
Rrp6, and hence with NNS. Consistently, the ncRNAs 
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most affected by loss of Setl contain more Nrdl/Nab3 
binding motifs and show higher Nrdl binding 
(Supplementary Figure S7D). Interestingly, NNS is 
recruited through interaction with Ser5 phosphorylated 
PolII C-terminal domain (7,8,10), and loss of Setl was 
recently shown to increase Ser5 phosphorylated 
C-terminal domain as well as Nab3 and Senl recruitment 
(49). Together, these observations suggest a model in 
which Setl and/or H3K4 trimethylation contribute to 
antisense production by interfering with NNS-mediated 
transcription attenuation. We speculate that this positive 
effect of Setl on asRNA transcription may be preferen- 
tially revealed in Arrp6, as our data indicate that the 
absence of Rrp6 compromises NNS early termination 
(31) and (Figure 4C). 

Promoter repression by Rrp6- and Nrdl-sensitive asRNAs 

Our study expands previous findings on the links between 
histone modifications and antisense production, and 
suggests potential mechanisms of how asRNAs contribute 
to silencing. Antisense-induced repression of the 28 Class 
(I) genes is strongly dependent on the presence of Setl, 
Hdal and Rpd3 (Figures 1 and 2). Setl -dependent H3K4 
dimethylation of non-coding upstream or antisense tran- 
scription has previously been implicated in gene repression 
by acting as a signal to recruit the Rpd3 or Set3 HDACs 
(24-26,28,50-52). Thus, H3K4 tri- and dimethylation de- 
posited by Setl during non-coding transcription may con- 
tribute to silencing in two ways, respectively, through 
stimulation of antisense production and by signalling the 
recruitment of specific HDACs to the sense promoter. 

Our analysis of promoter structure in the context of 
antisense classes suggests a link between asRNA produc- 
tion and promoter remodelling. The promoters of Class (I) 
and (IV) genes, which produce asRNAs and show dere- 
pression in the absence of Setl, Hda2 and Rpd3, are 
enriched for TATA boxes (as well as closed promoter con- 
figurations in the case of Class I). These genes are often 
transcribed in bursts, and the promoters tend to be 
regulated through extensive nucleosome rearrangements 
on the on-off switch (36,53), consistent with their sensi- 
tivity to histone modifications. Of note, our recent studies 
on the PH084 Class (I) gene using single-molecule FISH 
indicate anti-correlation of sense and asRNA expression 
in individual cells (31). Together these observations 
suggest a model in which, once a transcription burst is 
over, low-frequency antisense transcription through the 
promoters of Class (I) and (IV) genes contributes to 
nucleosome repositioning and occlusion of TF binding 
sites, thereby raising the threshold of the transcription 
on-switch (53). The HDACs Hdal and Rpd3 may either 
facilitate nucleosome repositioning or act downstream to 
consolidate repression by further reducing promoter 
accessibility. 

We showed that Nrdl depletion converts a 'non-func- 
tional' into a 'functional' asRNA (Class III), which is able 
to repress sense mRNA expression. The data suggest that 
the repression is linked to the frequency of asRNA mol- 
ecules that escape early termination and nuclear degrad- 
ation, so that they can extend into the sense promoter 



region. The transcription process per se over the 
promoter could trigger chromatin changes altering nucleo- 
some occupancy and/or histone modifications. The almost 
complete rescue of sense repression observed in Hda2 
mutants, even in the presence of high levels of antisense 
transcription over the promoter in the Nrdl-AA strain, 
highlights the crucial importance of histone deacetylation 
in the establishment or consolidation of the silencing 
process. Nearly 15% of the genes that accumulate anti- 
sense NUTs following Nrdl depletion show mRNA re- 
pression (14). It will be interesting to define what 
fraction of these responsive genes is repressed in a 
process dependent on Hda2 or through alternate mechan- 
isms based on transcription interference or polymerase 
collision (54). 

Regulation of antisense-mediated gene repression 

Our results indicate that a significant number of antisense 
transcripts are regulated at the level of early termination 
through a process implicating NNS and Setl. The recently 
published data from the Cramer laboratory confirm the 
global effect of Nrdl early termination on restricting 
ncRNA transcription and potential negative effects on 
gene expression (14). Our results indicate that the effi- 
ciency or strength of asRNA early termination may be 
variable, raising the possibility that leaky early termin- 
ation represents a selective advantage in the case of 
Class (I) and (II) genes, possibly because it may contribute 
to optimal regulation of gene expression. Interestingly, 
recent studies show that NNS activity is downregulated 
in response to changes in physiological conditions (55), 
potentially influencing asRNA expression levels and 
hence asRNA-dependent silencing. Notably, Rrp6 also 
undergoes physiological control and is downregulated 
during the onset of the meiotic program, giving rise to 
specific ncRNAs, that may be required for the proper ex- 
pression of genes at the entry of meiosis (56). Our obser- 
vations predict that these meiosis-specific ncRNAs result 
from increased read-through transcription of ncRNAs 
with weak early termination signals. 

A recent study indicates that co-transcriptional 
decapping of nuclear asRNAs by Dcp2 followed by 5—3' 
degradation by the nuclear exonuclease Rati may repre- 
sent an alternative mechanism to regulate antisense tran- 
scription and expression of the corresponding gene (27). 
We found little overlap between the 100 identified Dcp2- 
and the Rrp6-sensitive asRNAs described in our study, as 
most correspond to the unclassified genes (data not 
shown), suggesting that asRNA production may be 
regulated by at least two pathways depending on the 
gene: either nuclear decapping and degradation by Rati 
or early termination by NNS linked to degradation by the 
TRAMP-exosome complex. A fraction of Dcp2- and 
Rrp6-sensitive asRNAs have also been described as 
XUTs, non-coding transcripts that accumulate in the 
absence of the cytoplasmic 5—3' exonuclease Xrnl, 
indicating that they are exported and degraded in the cyto- 
plasm (27,31,52,57). It will be interesting to determine the 
nucleocytoplasmic distribution of asRNAs from different 
gene classes in different mutant backgrounds by 
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visualization of single transcripts, as we recently did for 
PH084 (31). Such experiments may reveal the relationship 
between asRNA transcription regulation, localization and 
effect on gene expression. The conservation of sense/anti- 
sense pairs through evolution and their wide occurrence in 
higher eukaryotic species (17,48,58) suggest that similar 
mechanisms may underlie the regulation of antisense pro- 
duction and their effects on the fine-tuning of gene expres- 
sion in metazoans. 
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